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We evaluate the cross section for pp — > £ + £~ir° in the forward direction and for large lepton 
pair invariant mass. In this kinematical region, the leading-twist amplitude factorises into a short- 
distance matrix element, long-distance dominated antiproton Distribution Amplitudes and proton 
to pion Transition Distribution Amplitudes (TDA) . Using a modelling inspired from the chiral limit 
for these TDAs, we obtain a first estimate of this cross section, thus demonstrating that this process 
can be measured at GSI-FAIR. 
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Transition Distribution Amplitudes (TDAs) [T] are 
universal non-perturbative objects describing the transi- 
tions between two different particles ( e.g. p — ► tt, tt — ► 7, 
7r — * p) . They appear in the study of backward electro- 
production of a pion [2], of 7*7 — * pir and 7*7 — * tttt 
reactions [5] as well as in hard exclusive production of a 
7*7r pair in pp annihilation: 

p{Pp)p{p P ) -> l*(q)ir(pn) -> i + (Pe+)Z~ (Pe-)^(p^) (1) 

at small t — (p v — p p ) 2 (or at small u = (p n — Pp) 2 ), 
which is the purpose of the present work. The TDAs are 
an extension of the concept of Generalised Parton Distri- 
butions (GPDs), as already advocated in [4]. The proton 
to meson TDAs are defined from the Fourier transform of 
a matrix clement of a three-quark-light-cone operator be- 
tween a proton and a meson state. They obey QCD evo- 
lution equations which follow from the renormalisation- 
group equation of the three-quark operator. Their Q 2 
dependence is thus completely under control. 

Whereas in the pion to photon case, models used for 
GPDs [13 13 13 H] could be applied to TDAs since they 
are defined from matrix elements of the same quark- 
antiquark operators, the situation is clearly different for 
the nucleon to meson TDAs. Before estimates based on 
models such as the meson-cloud model [3] become avail- 
able, it is important to use as much model-independent 
information as possible. In j2], we derived constraints 
from the chiral limit on the TDAs p — > n and made a first 
evaluation of the cross section for the backward electro- 
production of a pion in the large-£ (or small E v ) region. 
Related processes were also recently studied in [10] sim- 
ilarly to what was proposed in |llj . In this work, we 
apply the same setting to evaluate the cross sections for 
pp — > £ + £~ir° in the kinematical region accessible by GSI- 
FAIR [12 in the forward limit and at moderate energy 
of the meson. 

In the scaling regime where Q 2 = q 2 is of the or- 
der of W 2 — (pp + p P ) 2 , the amplitude for the pro- 
cess at small t - or CM angle of the pion 
close to - involves the p -> tt TDAs V p7r (xi, £ , A 2 ), 



A pv (xi,F,A 2 ), T p7T (xi,^,A 2 ), where x, (i = 1,2,3) de- 
note the light-cone-momentum fractions carried by par- 
ticipant quarks and £ is the skewedness parameter such 
that 2£ = x\ + X2 + X3. The amplitude is a convolution 
of the antiproton DAs, a perturbatively-calculable-hard- 
scattering amplitude and the p — > n TDAs. 




FIG. 1: The factorisation of the annihilation process pp — > 
7*7r into antiproton-distribution amplitudes (DA), the hard- 
subprocess amplitude (Mh) and proton — » pion transition dis- 
tribution amplitudes (TDA) . 

The momenta of the subprocess pp — > 7*7r are defined 
as shown in Fig. [T] The z-axis is chosen along the collid- 
ing proton and antiproton and the x—z plane is identified 
with the collision or hadronic plane. We define the light- 
cone vectors p and n such that 2 p.n = 1, as well as 
P = [p p + Ptt)/2, A — p v — p p and its transverse com- 
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ponent A T (A^ < 0). £ is defined as £ 
express the particle momenta through a Sudakov decom- 
position : 
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where a = IF 2 - 2M 2 + Wy/W 2 - AM 2 ~ 2PF 2 ; the 
approximate expression for g is obtained with M -C W. 
For At = 0, M <C W and m w = 0, one gets 

W 2 



p P = (i + Op 

2£M 2 



f = 



pp 



l 



Q 2 



K = (i - Op- 



(4) 



2 IF 2 - Q 2 



In the fixed-target mode, the maximal reachable value 
for W 2 = 2M 2 + 2ME p at GSI will be ~ 30 GeV 2 (for 
Ep = 15 GeV). The highest invariant mass of the photon 
could be Q 2 max ~ 30 GeV 2 . We refer to Ref Q3] for a 
complete discussion of the kinematically allowed domain. 
In terms of our notations, in the proton rest frame, we 
have p = j^(l,0,0,-l) and n = i±f (1,0,0, 1). Thus 
£ S [0.5, 1] corresponds to |p*| < M/3 ~ 310 MeV in the 
laboratory frame at At = 0. 

Let us now turn to the kinematics of pp — > ^ + ^~7r°. 
In general, we have for the unpolarised differential cross 
section: 



da = 



1 



2 v / A(W 2 ,M 2 ,M 2 )(2tt) 5 



|M| 2 d 3 (P,S). (5) 



The 3-particle differential Lorentz invariant phase space 
(dLIPS), ds(PS), can be decomposed into two 2-particle 
dLIPS (where q is the momentum of the 7*): 
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In the pp CM, we have: 

, d 3 ^ d 3 q 
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and in the £ + £ CM , we have (mt — 0): 

d 3 pi+ d 3 j>(- dQi dcos 
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Pp) 2 in terms of cos 9* [TJ, we get 
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VM^ 2 , M 2 : M 2 )^/\(W 2 ,Q 2 ,ml) (9) 



Altogether, by integrating on tp* and on tpi , 
da Jd^|7WPP^ + ^-°| 2 



dtdQ 2 dcos9 e MW 2 {W 2 - 4Af 2 )(27r) 4 
to be compared with the cross section for pp — > 7*71"° 
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dt 16ttW 2 (W 2 -4M 2 )' 

At A T = 0, the leading-twist TDAs for the 
p 7T° transition, V; p7r (x,, £, A 2 ), Af^, £, A 2 ) and 
Tf v (xi, £, A 2 ) are defined as (see Appendix for details) : 

^(<7r°( P7r )| e^u^mj^^njii^sn) |P(p P ,s p )) 

* JV yP 7r °(^C') Q/ 3(u + (p p ,Sp)) 7 

+ Ar°(^7 5 C)a /3 (7 5 w + (p P ,^)) 7 
+ T i P7r (c r PAl C) Q/3 (7' i M + (pp,Sp)) 7 

where cr^ v = 1/2[7 M , 7"], C is the charge conjugation 
matrix, f n = 131 MeV is the pion decay constant and 
/at ~ 5.2 • 10~ 3 GeV 2 . u + is the large component of 
the nucleon spinor : u(j> p ,s p ) = (jjip 1 + fof)u(p p , s p ) = 

u~(p p ,s p ) + u + (p p ,Sp) with u + (p p ,s p ) ~ Wpp and 



u (p P , s p ) ~ \Jl/p P 



For the three TDAs V pv , A*" 1, and T p7r , contributing 
in the limit Ay — > 0, we use the following expressions for 
Ay = and large £ (see Appendix): 



{Vr ,A? ,Tf }( Xl ,x 2 ,x 3 ,Z,A 2 ) 



(13) 



i{n^3^ } (|,|,|), 

where and T p are the proton DAs [TS]. 

At the leading order in a s and at At = 0, the ampli- 
tude M s x pSp for p(p p , Sp)p(p p , s p ) 7*(g, A)7r°(p w ) reads 



with 5! pS * 



54/ T g4 

5,,+/ 



(14) 
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(15) 



where [dx] = dxida;2(ia;3(5(2£ — J2k x k) an d [dy] = 
dxidx2d2;3(5(l — J2k Vk)] the coefficients i? Q (a — 1, 14) 
exactly correspond to T a in [5] after the replacement 
— it — > ie due to the presence of the 7* in the final in- 
stead of initial state. Even though the TDA formalism 
can be applied at any value of £ (or E v ), we have for 
now at our disposal estimates for the p — > n TDAs only 



3 



at large £. In the following, we shall therefore limit our- 
selves to the computation of the cross section for reaction 
(JlJ in this region. At large £, the ERBL regime (xj > 0) 
covers most of the integration domain. Therefore it is le- 
gitimate to approximate the cross section only from the 
ERBL contribution, i.e. when the integration range of 
the momentum fractions is restricted to [0, 2£]. 

The differential cross section for unpolarised protons 



and antiprotons is calculated as usual using Eq. (Ill 
from the averaged-squared amplitudes, 

l^l 2 = jE^-(>0*- 
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|.Moo| 2 vanishes at the leading-twist accuracy, as in the 

nucleon-form-factor case. The same is true for | 2 

and |.Mo+ 1 2 1 etc., since the x and y directions are not 
distinguishable when A\ is vanishing. We then define 
\Mt\ 2 = \M^\ 2 + \MZZ\ 2 . 

To compute X, we need to choose models for the DAs 
and the deduced TDAs. For the sake of coherence with 
experimental data, we shall use reasonable parametrisa- 
tions of CZ [15] and KS [IB], which are both based on an 
analysis of QCD sum rules. For CZ, they are 

V p {xi) = tp as [11.35(a; 2 -|-x^)-|-8.822; 2 -1.68a;3-2.94], 
A p { Xi ) = <p as [6.72{x\-x\% (17) 
T p ( Xl ) = <^ s [13.44(x 2 + x 2 J+4.62x 2 + 0.84x 3 - 3.78], 

and for KS (which we used in Fig. [2]) 

- 6.72x3 -5.04], 
a*)], (18) 
f 0.84x3 -7.56], 



V p (xi) = ^ as [17.64(x 2 + x%) + 22.68x; 
A p (xi) = ^ as [2.52(^-x 2 ) + 1.68(x 2 
T p (xi) = (/9 as [21.42(x 2 + .x 2 ) + 15.12x: 



and we evaluate our model TDAs from Eq. (13 1. This 
gives X ~ 1.28 • 10 5 for CZ and I ~ 2.15 • 10 5 for KS; this 
yields an induced uncertainty of order 3 for our estimates 
of the cross section. In the following, we use a s = 0.3 as 
suggested in [15]. 

Altogether, we have the following analytic results for 
the dominant ERBL contribution: 
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From this, we straightforwardly obtain whose W 2 
evolution is displayed on Fig. [2] (a) at A T = for the two 
extreme values of meson longitudinal momentum where 
one may trust the soft pion limit, corresponding to p% = 
or \p^\ = M/3 in the laboratory frame = 1 or 1/2). 
For the process ([T]), the averaged-squared amplitude is: 

s p ,sp : X, A' 



-^£ AA '-^(Xf s _)*, 

Q2 Q2^ s pSp / 



with £ AA ' = e 2 Tr(^-^(A)jV^*(A')). Integrating on the 
lepton azimuthal angle ipg, we have 

,27re 2 (l + cos 2 t ) 
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FIG. 2: (a) Differential cross section da/dt for pp —t j*it° as 
a function a W 2 for = (lower curve) and \p%\ = M/3. 
(b) Differential cross section da / (dtdQ 2 ) for pp — > £ + £~ir° as 
a function of Q 2 for various beam energies. 



from which we get, via Eq. (10 1 and integrating over 8g 
the differential cross section displayed in Fig. [2] (b) . 

Although the cross sections are evaluated at At = 0, 
we do not anticipate any dramatic Ay-dependence of the 
TDAs below a few hundred MeV, so that our estimates 
are likely to be valid in a not-too-narrow At region. To 
evaluate a magnitude of the integrated cross section we 
take as an example the kinematical region with W 2 = 10 
GeV 2 in Fig. 2b accompanied by the Q 2 window 7 GeV 2 
< Q 2 < 8 GeV 2 , which corresponds to £ w 1/2 or to the 
pion momentum of the order 310MeV. Integrating over 
this Q 2 -bin and in a t— bin corresponding to At < 500 
MeV leads then to a cross section around 100 femtobarns. 
Such a cross section is sizable and seems to be accessible 
to experimental setups such as PANDA with the designed 
value of the FAIR luminosity. 

The calculations done till now and which involve the 
proton — > 7T° TDA are valid for the small t region. Let 
us however stress that - due to the charge symmetry - 
an identical result will be obtained in the small u re- 
gion but with the p — * tt° TDA. In the laboratory frame 
at GSI-FAIR, this second region is quite different from 
the previous one since the 7r° meson is boosted in the 
forward direction. A precise detection of the particles 
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of the final state, either in the proton or the antiproton 
"fragmentation" kinematics, will depend on the detector 
performances in the respective regions. 

In conclusion, we have demonstrated that the study 
of proton-antiproton exclusive annihilation into a lepton 
pair and a pion is feasible at large values of the lepton 
pair invariant mass in the forthcoming PANDA or PAX 
experiments at GSI-FAIR. We believe that such a study 
will bring unique information about the inner structure 
of the proton, and particularly about the pion content 
of a proton, provided that the predictions (scaling be- 
haviour, angular dependence of the lepton pair) of the 
factorised framework used here are shown to be valid. 
Expressing the cross section of this process as the convo- 
lution of a hard-scattering amplitude for quarks with a 
photon with DAs and proton to pion TDAs will allow to 
evaluate these new hadronic matrix elements which con- 
tain much information about the confinement dynamics. 
It has to be emphasised that the same hadronic matrix 
elements appear also in a similar description of backward 
electroproduction of a pion [2] . 

Note that other channels are also of much interest, such 
as pp — ■> or pp — > l + l~ p° . The theoretical treat- 

ment of the j] case is identical to the one for the ir° case, 
but for the isosinglet nature of r\. In the p° case, one 
should distinguish between the longitudinally polarised 
meson where the TDA has the same structure as for the 
7T° and the transversally polarised case which leads to 
more TDAs. Needless to say, we are strongly lacking of 
model estimates for these p —> rj and p — > p TDAs but 
their experimental determination (or at least the mea- 
surement of their ratios to the p — > 7r TDAs ) opens a 
fascinating window on the properties of the sea quarks in 
the proton wave function. 
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is supported by the Polish grant 1 P03B 028 28, the 
French MAEE Eco-Net program, the EU contract FJI3- 
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Appendix. We now derive the general limit of the 
three contributing TDAs at At — in the soft-pion limit, 
when £ gets close to 1. In that limit, the soft-meson 
theorem [TT] derived from current algebra applies [TT] . 
which allows us to express these 3 TDAs in terms of the 
3 Distribution Amplitudes (DAs) of the corresponding 
baryon. Conventionally |15j . the three proton DAs are 
defined through the decomposition of the following ma- 
trix element of the 3-quark operator in terms of three 
invariant functions of the scalar product of the light- 
like separation Zin = Zi with the proton momentum p pi 
V p {S l .p p ),A p (S l .p p ) and T p (zi.p p ), 



(0\u a (z 1 n)u l3 (z2n)d J (z 3 n)\p p )) = |/ivX 
V p (z l .p p )(p p C) al3 (j 5 u + (p p , s p )) 7 



+A p (zi.p p )($ p -f b C) alj u + (p p , s p ) 1 
+T p (z l .p p )(a Pp ^ C) Q/3 (7 A1 7 5 w + (pj,, s p ))- 



(21) 



The latter functions satisfy V p (zi.p p = 0) = T p (zi.p p — 
0) = 1 and A p (zi.p p = 0) = 0, which provides the 
interpretation of /jy as the value of the proton wave 
function at the origin. To go to momentum space one 
writes a Fourier transform [16 which enables to define 
functions of momentum fractions Xi (F — V P ,A P ,T P ) 
([dz.-] = d(z 1 .-)d(z 2 .-)d(z 3 .-)) : 



F(xi) = 



[dz.pp] 
(2tt) 3 



e^^^-^F^i.pp). 



(22) 



Inspired by [TT] , which considered the related case of the 
distribution amplitude of the proton-meson system, we 
use the soft pion theorems [17] to write: 



(■K a {p«)\0\P{ PllSl )) 



(0|[Q»,O]|P( Pl , Sl )), (23) 



where we neglected the nuclcon pole term, which does 
not contribute at threshold. 

For the transition p — > 7r°, Qf = Q\ and O — u a upd T 
Since the commutator of the chiral charge Q5 with the 
quark field ip (r a being the Pauli matrix) 



-7 V> > 



(24) 



the first term in the rhs of Eq. ( |23| gives three terms 
from {^ b u) a Ufjd 1 , u a {^ b u)pd 1 and u a Uf3{^ b d) 1 . The cor- 
responding multiplication by 7 s (or (7 s ) T when it acts 
on the index (3) on the vector and axial- vector structures 
of the DA (Eq. (21 1) gives two terms which cancel and 



the third one, which remains, is the same as the one for 
the TDA up to the modification that on the DA decom- 
position, p p is the proton momentum, whereas for the 
TDA one, p is the light-cone projection of P, i.e. half 
the proton momentum if one neglects the pion one. This 
introduces a factor 2£ in the relations between the 2 DAs 

A p and V p and the 2 TDAs and Af °. 

To what concerns the tensorial structure multiplying 
T p , the three terms are identical at leading- twist accu- 
racy and correspond to the structure multiplying T p7T , 
this gives a factor 3. We eventually have the soft limit 
for our three TDAs at A T = 0: 

(vr°,Af,Tr a m. P ) = ^^sno^) ( 25 ) 

We will derive now this relation in the momentum rep- 
resentation of DAs and TDAs. To do so, we start from 
translational invariance which implies 



\u a {z\ + a)up(z 2 + a)d^(z 3 + a)\p p ) = 
e- ia - p "(0\u a (S 1 )u p (S 2 )d 7 (S 3 )\p p ), 



(26) 



and thus F((zi + a).p p ) = e la ' Pp F(zi.p p ). In momentum 
space, we correspondingly get 



F{ Xi ) = 



[d(z + a).p p ] lY . k x k (z k +a). Vv 



(27T) 



i(£ k x k -l)a.p p 



F( Xi ). 



F((zi + a).p p ) 
(27) 
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This condition is conveniently expressed by the following 
redefinition: F{xi) = 5(^2 k x k — l)F(xi). The inverse 
Fourier transform is then written as ([dx] — dxidx2dx 3 ): 



F{zi-P P ) 



[dx] e-^ kXkik - p ?5Cy x k - l)F( Xi ). (28) 



k 



The normalisation conditions then reads 



The condition derived from translational invariance 
would then be satisfied by S(J2k x k a -Po — a.(PA — Pb))- 
In order to get the correct support, i.e. <5(2£ — J2k x k)> 
for a translation a along n, we have to choose po such 
that 2£ = "- (pa ~ Pb) = which is satisfied by 

Po = (pa + pb)/2 = P or a light-cone vector p (p 2 = 0) 
such that P.n — p.n. We choose p and have 



f\dx] S(£ x k - 1)(V* A", T>)( Xi ) = (1, 0, 1). (29) m " £ X *) F *-+B(*i) = F{F{zi{n.p))) = (31) 

oo 



Note that the delta function insuring translational in- 
variance is exactly the one expected from the interpre- 
tation that Xi be the momentum fraction carried by the 
quark i off the proton of momentum p p . This shows 
that the natural conjugate variable to the momentum 
fractions the Zi.ppS, i.e. the spatial separation 

dotted by the proton momentum. Indeed, p p enters in 
the exponential of the rhs of Eq. (27 1, via the conjugate 



variable to Xi, Zi.p p , and as the initial-state momentum 
of the matrix element. Another choice than Zi.p p would 
not have provided the correct support for the x^s. 

The case of TDAs is similar except for the choice of the 
natural conjugate variable. We start with translational 
invariance: 

{B{ PB )\t ijk u l a {z l + a)u^(z 2 + a)d k (S 3 + a)\A(p A )) = 

e ---(^-P B ) (jB(pB) | e yv a (5 1 )^(^)^(z 3 )|A(^)). 

Now let us define a Fourier transform without specifying 
the momentum po to which we dot the spatial separation: 



F A ^ B (xi) 



' J [d((z + a).po)} 
o (2^) 3 



e iS k x k (z k +a).p x /-gg^j 



F A ^B((zi + a). Po ) = e l ^ Xka ^- a ^ A 'P^F A ^ B ( Xl ). 



{n.pf 



[dz] 
(2^)3 



3 iSkXkZ ^ n - p) F(zi(n.p)). 



Let us now use the Fourier transform Eq. ( 22 1 on both 



side of Eq. (251. Defining a such that p p .n = aP.n = 
ap.n, we get for instance for the V's: 



[dz.pp] 
(2tt) 3 



/ ^e^ kXkZka(n ^Vr° ( Zl (n.p)) 



= a 3 5(2t-J2(axkWr°(ax l ). 



(32) 



We conclude that translational invariance imposes natu- 
rally, through the delta function, that a — 2£; the change 
of variable x\ — axi then yields 



_0 



{vr ,a? ,rr }(x u x 2 ,x 3 ,z,a 2 ) 



(33) 
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